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Abstract 

The dynamics of nonlinear reaction-diffusion systems is dominated by the onset of patterns 
and Fisher equation is considered to be a prototype of such diffusive equations. Here we inves- 
tigate the integrability properties of a generalized Fisher equation in both (1+1) and (2+1) 
dimensions. A Painleve singularity structure analysis singles out a special case (m = 2) as 
integrable. More interestingly, a Backlund transformation is shown to give rise to a linearizing 
transformation for the integrable case. A Lie symmetry analysis again separates out the same 
m — 2 case as the integrable one and hence we report several physically interesting solutions 
via similarity reductions. Thus we give a group theoretical interpretation for the system under 
study. Explicit and numerical solutions for specific cases of nonintegrable systems are also 
given. In particular, the system is found to exhibit different types of travelling wave solutions 
and patterns, static structures and localized structures. Besides the Lie symmetry analysis, 
nonclassical and generalized conditional symmetry analysis are also carried out. 

1 Introduction 

Integrable systems play the role of prototypical examples to identify and understand the vari- 
ous phenomena underlying nonlinear dispersive systems such as Korteweg-de Vries, sine-Gordon, 
Heisenberg spin, nonlinear Schrodinger, Davey-Stewartson, Kadomtsev-Petviashivili, etc. equa- 
tions [Ablowitz & Clarkson, 1991; Lakshmanan & Rajasekar, 2003]. Nonintegrable perturbations 
can be analysed in terms of the basic excitations of such integrable models [Scott, 1999]. These 
integrable systems are often shown to be linearizable in the sense that they can be associated 
with two linear differential equations (of which one is a linear eigenvalue problem), namely the 
so called Lax pairs. In (1+1) dimensions, the exponentially localized stable entity, namely the 
soliton, turns out to be the basic structure exhibiting elastic collision property. The soliton excita- 
tion has remarkable physical and mathematical properties. In particular, the underlying nonlinear 
evolution equations exhibit infinite number of Lie-Backlund symmetries [Bluman & Kumei, 1989]. 
The Lie point symmetries lead to similarity reductions which are related to Painleve type ordinary 
differential equations (ODEs) which are free from movable critical singular points. More generally, 
the integrable soliton equations satisfy the Painleve property and the solutions are in general free 
from movable critical singular manifolds. 

In the case of nonlinear reaction diffusion equations, the dynamics is dominated by the onset 
of patterns [Murray, 1989; Walgraef, 1996; Lakshmanan & Rajasekar, 2003]. Out of all possible 
modes available to the system, they tend to select the most stable structures which give rise to 



various patterns. These patterns range from simple to chaotic, depending on the nature of diffu- 
sion, nonlinear reaction terms and external forces as well as the spatial dimension. Some of the 
dominant patterns are 

(1) homogeneous or uniform states 

(2) travelling waves 

(3) wavefronts and pulses 

(4) Turing structures: a) stripes b) spirals c) scrolls, etc. 

(5) spatiotemporal chaos 

and so on. Thus mode selection and stability dominates the study of such systems. Some of 
the well known nonlinear diffusion and reaction-diffusion systems [Murray, 1989; Walgraef, 1996; 
Whitham, 1974] include Burgers' equation, Fisher equation, Kuromoto-Sivashinsky equation, 
Gierer-Meinhardt equation, FitzHugh-Nagumo equation, Belousov-Zhabotinsky reaction equation, 
Brusselator model equation and so on. A large amount of literature on such systems is available. 

In view of the complexity involved in analysing nonlinear diffusion equations, it will be very 
valuable to identify integrable nonlinear diffusive systems and to build on them the study of 
nonintegrable cases. Burgers' equation is a standing example of such an integrable case [Sachdev, 
1987]. It is linearizable in the sense that the Burgers' equation 

u t = Du xx - uu x , (1) 

under the Cole-Hopf transformation 

u = -2D^- (2) 

v 

gets transformed into the linear heat equation 

v t = Dv xx . (3) 
Similarly the Fokas-Yortsos-Rosen equation [Fokas & Yortsos; 1982 & Rosen, 1982] 

u t = u 2 (Du xx - u x ) (4) 
under the variable transformations [Rosen, 1982] 

x = — Dhxv, u = — D— , (5) 

v 

gets again transformed into the linear heat equation Eq. Thus both the above systems are 
linearizable and they may be considered to be C-integrable in the Calogero [1991] sense. Also both 
the Burgers' equation and Fokas-Yortsos-Rosen equation possess interesting Lie point symmetry 
structures and infinite number of Lie-Backlund symmetries. So, it will be of interest to investigate 
whether other integrable and linearizable nonlinear reaction-diffusion equations exist and if so 
what is the role of symmetry therein and what kind of solution structures and patterns exist in 
them. 

Considering nonlinear reaction-diffusion equations, Fisher equation is considered to be proto- 
typical [Murray, 1989]. Typically, its form reads as 

u t - Am - u + u 2 = 0, (6) 

where A is the Laplacian operator. In one dimension it admits travelling wavefronts with velocities 
c > c m i n — 2, depending on the initial condition [Murray, 1989] and in two dimensions interesting 
spatiotemporal patterns arise; however it is not integrable in any dimensions. In one dimension, 
however it is known to have an exact propagating wavefront solution of the form [Ablowitz & 
Zeppetella, 1979] 

-2 

(7) 

where k is an arbitrary constant, corresponding to the wave speed c = v5/6. No other exact 
solutions are known for Eq. ©. 



u(x, t) = 1 — 



1 + 7i exp 



V6 
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2 Generalized Fisher (GF) Equation 



In this connection, it has been pointed out [Grimson & Barker, 1994] that there arises interest- 
ing generalization of the Fisher equation in the description of bacterial growth in multicolony 
environment, in chemical kinetics and in various living phenomena. Its general form is 

= DAu(r,t) + A{u)[Vu(r,t)) 2 + Xu{r,t)G{u), (8) 

where D is the diffusion coefficient, V and A are gradient and Laplacian operators respectively. 
To be specific, Eq. JSJ represents the spatiotemporal growth of bacterial colonies with both local 
and nonlocal modes of growth. The third term on the right-hand side of Eq. © represents the 
local growth characterised by the local growth rate A and a local growth function G(u). Suppose 
the colony is unable to add cells locally due to packing constraints, the nonlocal addition of cells 
on the surface of the colony occurs to account for the amount of growth required by the global 
growth law. Such an addition occurs as a consequence of the expansion and reorganization of the 
colony and is represented by the second term while the first term represents diffusion. 
An important special case of Eq. (JHJ) is the generalized Fisher type equation 

u t - Au- j^— (Vu) 2 - u{\ - u) = 0, (9) 

where the subscript denotes partial differentiation with respect to time. In the study of population 
dynamics, u(r, t) refers to the population density at the point r at time t. In Eq. (|5J), the nondiffu- 
sive linear term modelling the birth rate gives rise to an exponential growth of it(r, i) in time while 
the quadratic term that models competition between individuals for food, etc. leads to a stable, 
homogeneous value u = 1 at long times and the diffusion term represents the spatial variation of 
the population. Further the third term represents the nonlocal growth of cells while the fourth 
term corresponds to the local growth. This introduces the possibility of spatial pattern formation 
between the homogeneous regions with u = 1 and u = for appropriate initial conditions. Note 
that for m — 0, Eq. 10 reduces to the standard Fisher equation ijfJJl. Again Eq. © is in general 
nonintegrable, except for the special case m = 2. It is linearizable under the transformation 

" = 1-^, (10) 

as pointed out by Wang et al. [1996] for the one dimensional case (see also Brazhnik and Tyson 
[1999a,b]). Under this transformation the generalized Fisher equation can be shown to reduce 
to the form 



Xt -A X -x = 0, (11) 

which is nothing but the linear heat equation in arbitrary dimensions. Thus every solution of l|llfl 
corresponds to a solution of Eq. via the transformation (|10l) subject to boundary and initial 
conditions. Note that the transformation holds good irrespective of the spatial dimensionality 
of the generalized Fisher system © as long as the parameter takes the value m = 2. Such a 
linearization in any dimensions is a rare situation indeed. Even in soliton theory, straightforward 
extensions of integrable equations lose their integrability as soon as the dimension is increased. In 
fact Wang et al. [1996] and Brazhnik and Tyson [1999a, b] make use of this linearizability property 
to obtain several interesting wave solutions, spatiotemporal patterns and static structures. 

The very special nature of the m = 2 case prompts us to investigate the integrability nature 
of the generalized Fisher equation @ from different points of view, particularly from singularity 
structure and symmetry points of view in order to understand why the case m = 2 alone is 
integrable. For the nonintegrable case (m ^ 2) we investigate special solutions. 

In Sec. we show that the PDE © in one dimension is free from movable critical singular 
manifolds only for the special case m = 2. Further we point out that a Backlund transformation 
gives rise to the linearizing transformation in a natural way for the integrable case. In Sec. 
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21 we derive certain explicit and numerical solutions for the GF equation in one dimension via 
symmetry analysis and similarity reductions. In Sec. we present the associated Lie algebra for 
the (2+1) dimensional GF equation. Then we report different types of propagating structures, 
static solutions and localized structures through Lie symmetry analysis in Sees. and {7\ In 
addition, exact solutions for special cases of the nonintegrable GF equation are also given. We 
give the nonclassical and generalized conditional reductions in Sees. |H1 and El respectively. Finally, 
we summarise our results in Sec. 1101 

3 Singularity Structure Property of the GF Equation 

It has been argued for sometime that the Painleve property (P-property) as proposed by Weiss 
et al. [1983] is a working test to identify whether a given nonlinear partial differential equation 
(PDE) is integrable or not. The Painleve property demands that the solution of the given non- 
linear PDE is free from movable critical singular manifolds (branching types, both algebraic and 
logarithmic as well as essential singular), so that the solution is single valued in the neighbourhood 
of a noncharacteristic movable singular manifold, </>(r,t) = 0, <fi x ,4> y ,(j)t ^ 0. The Weiss Tabor 
Carnevale (WTC) procedure then demands that the solution u(r, t) for Eq. © can be expanded 
around the noncharacteristic movable singular manifold as a Laurent expansion, 

oo 

u = Y j n ] {^tW +P (12) 

3=0 

where Uj's are functions of r and t and p is the leading order exponent of the Laurent expansion. 
For single valucdncss of the solution, one requires that p is an integer and that there exist sufficient 
number of arbitrary functions, Uj(r,t),j = 1,2, ...,N (for Eq. © it is two), so that the Laurent 
expansion H12|) can be considered as the general solution of the nonlinear PDE @ without the 
introduction of movable critical singular manifolds. Following the standard algorithm we can 
proceed as follows. As an example, we demonstrate the analysis for the one-dimensional case, 
which can be then extended to the higher dimensional case without difficulty. 

3.1 P-property of eq. |fP|) 

In the one-dimensional case, the generalized Fisher equation can be written as 

m 2 2 

u t - u xx - - u x — u + u = 0. (13) 

1 - u 

The WTC algorithm then identifies the following results. Starting with the leading order behaviour 
u sa u c/> p as suggested by the form (|12f> . one can identify three possibilities: 

(i) V = -2 

1 

(it) p = , m^l 

1 — m 

(in) p = 0. (14) 

One can easily check that for all the above three leading orders, only for the value m — 2, the 
solution is free from movable critical singular manifolds. Actually one finds that for m = 2 and 
p = — 1, the leading order coefficient uq is an arbitrary function in addition to the arbitrary nature 
of the singular manifold so that the Laurent series contains the required number of arbitrary 
functions in terms of which all the other coefficient functions Ui(x,t),i > 1, can be uniquely 
found. For the other two leading orders with m = 2, one identifies only one arbitrary function but 
without the introduction of any movable critical singular manifolds and so they may be considered 
as corresponding to special solutions. The above analysis also holds good in higher dimensions of 
@ as well. 
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Consequently one finds that the P-property is satisfied only for the special choice of parameter, 
m — 2 for Eq. @ and for all other values of 2) it fails to satisfy the P-property and the 
system is expected to be nonintegrable. Note however that there may exist certain special singular 
manifolds <fit — 4>xx = 0, for which the P-property may be satisfied in special situations. However 
in general for the special choice of the parameter m = 2 of Eq. © the P-property is satisfied for 
general singular manifolds and the system is expected to be integrable in this case. 

3.2 Backlund transformation 

Now we consider the Laurent series Ijl2(l applicable to Eq. I|13|l and truncate it at constant level 
term, by making all other coefficient functions Uj(r, t) to vanish for i > 1 consistently: 

u = — + iti. (15) 

We demand that if u\ is a solution of (|13() . so also u is a solution. Then, from Eq. Ijl5|) . we have 
u u - u lxx -ui(l- ui) + (u t - Uoxx ~ u {l - 2wi))0~ 1 + {2u 0x (j) x 



-u (4> t - 4>xx - u )4> - 2u <j> x 4> - 2(f) 



(1 - ui)4> - it 



+ 2u 0x u lx (t) 1 + (ul x - 2u u lx (j) x )(j) 2 - 2u UQ X (j) x <$> 3 + ulfy 4 =0. (16) 

Consequently one finds that Uq and 4> satisfy a set of coupled PDEs arising from Eq. (|16|) . 
Starting from the trivial solution, u\ = 0, of Eq. (fO|l . we find that the choice 

u = (b (17) 

solves Eq. H16|) . As a result, we find the new solution 

u=^+u x = l (18) 

which is indeed a solution to Eq. (|13fl . 

Taking the solution (|18|) as ui now, we find that the system l|16|) satisfies the equations 

m = -1, 

4>t -<l>xx-<t>+l = 0. (19) 
Now defining the new function \ as <fi — 1 + X one obtains the linear heat equation 

Xt ~ Xxx ~ X = 0. (20) 
Thus the transformation [Bindu et ai, 2001; Bindu & Lakshmanan, 2002] 

u = l--±—, (21) 
1 + X 

is the linearizing transformation for Eq. I|13|l , where the function \ satisfies the linear heat equation 
(|20|l . We note that this is exactly the transformation given by Wang et al. [1996] in an adhoc 
way. Here the transformation is given an interpretation in terms of the Backlund transformation. 
Thus the application of Painlcvc singularity structure analysis selects not only the integrable case 
associated with Eq. l|13fl . but also helps to identify the associated linearizing transformation. 

One can generalize the above analysis to higher dimensional case of the Eq. © for the choice 
m = 2 and arrive at the same general conclusion that the transformation (|21|l can be interpreted 
as a Backlund transformation. 
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4 Lie Point Symmetries of the (1+1) Dimensional GF Equa- 
tion and Integrability 

In this Sec, we investigate group invariance properties of the generalized Fisher equation in its 
(1+1) dimensional form given in Eq. (fO|l and discuss different underlying patterns via similarity 
reductions. First we consider the invariance of Eq. Ill.'iH under the one-parameter Lie group of 
infinitesimal transformations and derive the similarity reduced ODEs for integrable (to = 2) and 
nonintegrable (m ^ 2) cases separately. We solve the ODEs and provide explicit solutions. We 
show that the earlier reported solutions are sub-cases of our results. We construct the linearizing 
transformation from the infinitesimal symmetries and provide a group theoretical interpretation 
for it. 

4.1 Lie symmetries 

An invariance analysis of Eq. (|13[l under the one-parameter Lie group of infinitesimal transforma- 
tions 



t — > T = t + £r(i, x, u), x — ► X = x + e£(t, X, u), 

u — ► U = u + e4>(t,x,u), (22) 

singles out a special value of to, namely to = 2, to have nontrivial infinite-dimensional Lie algebra 
of symmetries 

r = a, £ = b, 4> = c(t,x)(l - uf, (23) 

where a, b are arbitrary constants and c(t, x) is any solution of the linear heat equation Ct — c xx — c = 
0. The infinitesimal symmetries <|23ll are obtained by following the usual procedure of solving the 
determining equations for the infinitesimal symmetries [Bluman & Kumci, 1989]. Throughout this 
paper we use the computer program MUMATH [Head, 1993] to determine infinitesimal symmetries 
and symmetry algebra. 

The Lie algebra of infinitesimal symmetries of the Fisher equation is spanned by two commuting 
vector fields and a vector field associated with an infinite-dimensional subalgebra, 

X 1= d t , X 2 = 8 X , X c = c(t,x)(l-u) 2 d Ul (24) 

respectively. In the above, the vector fields X\ and Xi reflect the time and spatial invariance of 
Eq. Ijl3|) . The commutation relations between these vector fields are given by 

[X!,Xa] = 0, [X 1 ,X c ] = X Ct , [A 2 ,A c ]=A Cx , (25) 

where X Ct = q(1 — u) 2 d u and X Ct = c x (l — u) 2 d u . For all other values of the parameter to(^ 2), 
Eq. H13|) admits only the trivial translational symmetries 

T = a, £ = b, = 0, (26) 

with the corresponding symmetry algebra 

[X 1 ,X 2 ]=0, (27) 

where X\ = dt and X 2 = d x . 

4.2 Similarity reductions 

We solve the characteristic equation associated with the symmetries (|23|) and (|26|l and obtain 
similarity variables in terms of which Eq. 1)13(1 can be reduced to an ODE. From the general 
solution of the ODE we construct explicit solutions for the GF equation fT3|) . 



G 



4.2(A) Integrable case (m = 2) 

Solving the characteristic equation 

dt 
a 



dx 

T 



du 



c(t,x)(l - u) 2 

associated with the symmetries (123(1 . we obtain the similarity reduced variables 

a 

z = ax — bt, u = 1 ;-r j— ■. — . 

a + v(z) + J c(t, x)dt 

In Eq. H29|) we can consider the quantity 

~(v(z) + / c(t,x)dt S j = x, 



(28) 



(29) 



(30) 



and verify that x satisfies the linear heat equation l|2U|l . Then the similarity transformation given 
by Eq. (|29|l is the linearizing transformation thereby giving a group theoretical interpretation of 
it. 

Making use of , Eq. {T3J| reduces to an ODE 



a 2 v" + bv' + v = 0, 



d 

dz ' 



whose general solution can be written as 



v = he m ^+I 2l 



m 2 z 



mi 2 = 



-b± V6 2 - 4a 2 
2a? ' 



(31) 



(32) 



Here I\ and 1% are integration constants. Now with the use of Eq. 129(1 . the solution to the original 
PDE (O becomes 



4a 2 > 0; 



4a' 



a + Iie mi(ax-bt) + he m 2 {ax - bt) + r c ^ x ) dt ' 
a 

a + eA ax - bt ) (Ji + I 2 (ax - bt)) + J c(t, x)dt ' 

a 

a + e A ax ~ ht ) (h cos q(ax - bt) + I 2 sin q(ax - bt)) + J c(t, x)dt ' 

b 2 - 4a 2 < 0, 



(33) 



with p = -b/2a 2 , q = V4a 2 - 6 2 /2a 2 . 

We note the following: By fixing the integration constants 1\ and I 2 suitably and for the choice 
c(t, x) — we can deduce all the interesting travelling wave solutions and stationary structures 
discussed by Brazhnik and Tyson [1999a, b]. For example, let us consider b 2 — 4a 2 > 0. Taking 
either one of the constants I\ or I 2 equal to zero we have a simple travelling wave solution. The 
choice I\ = 1% 7^ leads to a V-wave pattern. A Y-wave solution is produced by I\ = —I 2 . An 
oscillating wave front is obtained by restricting I 2 = with b 2 — 4a 2 < 0. The choice b 2 — 4a 2 = 
leads to an inhomogeneous solution when Ii = 0, 1 2 ^ 0. 

4.2(B) Nonintegrable case (m 7^ 2) 

A similar analysis on Eq. 1(26(1 for the (m ^ 2) case leads to the similarity variables 

z — ax — bt, u — w(z), (34) 
which reduces the original PDE 1(13(1 to an ODE of the form 

a 2 vv" — ma 2 v' 2 + bvv' — v 2 (l — v) = 0, v = 1 — w. (35) 
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For b ^ 0, this equation is in general nonintegrable except for m = and b/a — 5/y6- This special 
choice leads to the cline solution 



u{x, t) = 1 



1 H = exp 

\/6 



k ( x Ve* 



ia/6 



(36) 



where fc is an arbitrary constant reported by Ablowitz and Zeppetella [1979] (Fig. 1) and the 
surface plot is shown in Fig. 2. 



O 



1 O 



2 O 



Figure 1: An exact wavefront solution of the Fisher equation (£ = x — -jgt) 




Figure 2: Surface plot of the solution Ij36|l ■ 
Case (a): Static solutions (6 = 0, ra ^ 0) 

This is the time-independent (static) case and the similarity variables z = x, w = u lead to 
the reduced ODE of the form 

w" + -^-u/ 2 + tcfl - to) = 0. (37) 
1 — w 

This readily gives the first integral [Mathews & Lakshmanan, 1974] 

v' 2 = hv 2m + -^—v 2 - - 2 n v 3 , v = l-w, m^l, 3/2, (38) 
1 — m 3 — 2m 

where I\ is the integration constant. Eq. (|38|) on integration leads to elliptic or hyperelliptic 
function solutions for suitable values of m. For example, let us consider m = 1/2. Then Eq. 
becomes 

v' 2 = [h + 2v- v 2 ] v. (39) 
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Integration of Eq. I|39|l leads to 



dv 



y / v{I 1 +2v-v 2 ) 



(40) 



where I2 is the second integration constant. Let C\ 
Eq. l |%0*jl becomes 

dv 



^/(ci - v)(c 2 - v)(-v) 



1 + Vl + h, c 2 = 1 - VT+7T. Then 



(41) 



Now the LHS of the above equation can be integrated in terms of elliptic functions and the 
solutions are tabulated in Table 1. In a similar way one can integrate Eq. I|38[) and obtain elliptic 
and hyperelliptic function solutions for certain values of m. Besides the elliptic function solutions, 
for a particular value I\ = in Eq. i|39|) . we obtain a static solitary wave solution (Fig. 3) 



u = 1 - 



(3 - 2m) 
(2 - 2m) 



sech h — 



1 — m 



m < 1. 



(42) 



Eq. I|42|) is nothing but a limiting case of the elliptic function solutions. 



U 



-20 



2 



Figure 3: A static solitary wave pulse for m = 1/2 of the generalized Fisher equation, Eq. (|42|l . 
Case (b): Qualitative analysis 

In general Eq. I|35(l is of nonintegrable type and this leads to the use of numerical techniques 
to explore the underlying dynamics. Here we determine the equilibrium points and then study 
the system dynamics in the vicinity of these equilibrium points. Eq. (|35|l . after a rescaling, can 
be written as 



P 



P 



m Bp + v(l — v), ' = d/dz, 

v 



(43) 



where z = z/a and B = b/a. Then the equilibrium points are found to be (0,0) and (1,0). 
While the origin is an unstable equilibrium point, the stability determining eigenvalues for the 
equilibrium point (1,0) are found to be 



Ai. 



-B ± ^B 2 - 4 



(44) 



Now the nature of the singular point (1,0) is investigated as a function of the control parameter 
B in the range (—00, 00) by analysing the form of the eigenvalues A^a given by Eq. (|44p> and the 
results are tabulated in Table 2 and the corresponding plots are shown in Fig. (4). One identifies 
the existence of periodic pulses and wave front solutions in this case. 
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Figure 4: Propagating patterns and corresponding phase portraits in the v — v' plane of Eq. 13511 
(a) Periodic pulses, (b) Center, (c) Travelling pulse, (d) Stable focus, (e) Travelling wavefront, (f) 
Stable node. 

Table 1: Table 2. Classification of equilibrium point (1,0) of Eq. (|43[) as a function of the 
parameter B 



S.No. 


Range of B 


Nature of eigenvalues 


Type of attractor/repeller 


1. 


-co < B < -2 


Ai, 2 > 0, Ai ^ A 2 


unstable node 


2. 


B = -2 


A1.2 > 0, Ai = A2 


unstable star 


3. 


-2 < B < 


A1.2 = a Hj3, a > 


unstable focus 


4. 


B = 


Ai, 2 = ±i, 


center 


5. 


< B < 2 


Ai,2 = a ± a < 


stable focus 


6. 


B = 2 


Ai,2 < 0, Ai = A 2 


stable star 


7. 


2 < £ < oo 


Ai, 2 < 0, Ai ^ A 2 


stable node 



4.3 Relation between symmetries of m = 2 case and of linear heat 
equation 

In subsection 13.21 we derived a Backlund transformation that transforms the GF equation to a 
linear heat equation via singularity structure analysis. In this Sec, by employing the ideas given 
by Bluman and Kumei, [1989] (Theorem 6.4.1-1,2, p. 320), we derive the linearizing transforma- 
tion from the symmetries and thereby confirming the group theoretical interpretation for it. The 
aforementioned theorem gives a necessary and sufficient condition for the existence of an invertible 
mapping which linearizes a system admitting an infinite-parameter Lie group of transformation. 
That is, if a given nonlinear system of PDEs admit an infinite-parameter Lie group of transfor- 
mation having an infinitesimal operator 

X = Zi(x,u)d Xi +if(x,u)a u », (45) 
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with 

m 

£i(z,u) = S ^a1(x 1 u)F a {x 1 u), (46a) 

<T=1 

tn 

rf{x,u) = Y,fi(x,u)F"(x,u), (46b) 



(7=1 

then there exists a mapping 

= <f>j(x,u), j = 1,2, ...,n (47a) 
w 7 = V 7 (^,w), 7 = l,2,...,m (47b) 

that transforms the given system to a linear system of PDEs. Here F a is any arbitrary solution 
to the linear system of PDEs. In such a case, the components of the mapping cf>j and ip 1 satisfy 
the follwing set of PDEs: 

al{x,u)— +(3Z(x,u)— = 0, ir=l,2,...,m (48a) 
— +/^( a; ,t t )_ = <5 7 °", cr=l,2,...,m (48b) 

where <5 7CT is a kronecker symbol, 7, a = 1, 2, to. 

In our case, the infinitesimal operator generating the infinite parameter group of transformation 

is 

X = X c = c(t,x)(l-u) 2 ^-. (49) 

From Eqs. l(l8")l and l(4^|) we identify 

aj = a i = 0, ft 1 = (1 - u) 2 . (50) 

Then Eq. |@HK ) for <& becomes 

$„ - 0. (51) 
Clearly two independent solutions of Eq. (|51l) can be chosen as the new independent variables 

z\ =t, z 2 = x. (52) 

Equation H48b~) for becomes 

(l-u) 2 V u = l. (53) 
A particular independent solution of Eq. i|53|l can be easily found to be 

1 

w = . 

1 — u 

As a result we obtain a mapping 



(54) 



Zi=t, z 2 — x, w = (55) 

1 — u 

which transforms the Fisher equation to a linear heat equation 

v Zl - v Z2Z2 —v = 0, w = 1 + v. (56) 
Here we note that the above transformation is the same as the Backlund transformation (|21|l . 
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5 The (2+1) Dimensional GF Equation 



In general models that admit exact solutions are of considerable importance for understanding 
general behaviour of nonlinear dissipative systems. In one dimensions such models have received 
considerable importance. However, many realistic models are of two or three dimensional in nature 
and in this direction Brazhnik and Tyson [1999a,b] considered Eq. J^J in two spatial dimensions 
and explored five kinds of travelling wave patterns, namely, plane, V- and Y-waves, a separatrix 
and a space oscillating propagating structures. Thus, it will be quite interesting to study the role 
of symmetries that allows the system to exhibit different spatiotemporal patterns and structures 
which usually will possess some kind of symmetry. Moreover, the GF type equation does not 
lose its integrable nature (for to = 2 case) on its dimension being increased. This is a very rare 
situation and so we consider the generalized Fisher equation in 2-spatial dimensions, 

u t - u xx - u yy - - — (ul + Uy) - u + u 2 = 0, (57) 

and explore the underlying patterns. 

Extending the Lie symmetry analysis to Eq. (|57J) we find that the infinitesimal transformations 

t — ► T = t + sr(t, x, y, u), x — ► X = x + e£(t, x, y, u), 

y — > Y = y + £T](t, x, y, u), u — > U = u + e<j){t,x,y,u), e <C 1, (58) 

again singles out the integrable case for m = 2 with the symmetries 

t = a, £ = by + l, i] = —bx + n, <f> = c(t, x, y)(l — u) 2 , (59) 

where c(t, x, y) is any solution to the two dimensional linear heat equation ct — c xx — c yy — c = 
and b, I and n are arbitrary constants. For other values of 2), we get the following symmetries 

t = a, £ = by + I, rj — —bx + n, = 0. (60) 
5.1 Lie algebras 

The symmetry algebra of the (2+1) dimensional Fisher type equation Ij57(l with m = 2 is spanned 
by four vector fields 

Xi = d t , X 2 =yd x -xd y , X 3 =d x , X 4 = d y , (61) 

and a vector field associated with the infinite-dimensional subalgebra 

X c = c(t,x iy )(l-u) 2 d u . (62) 

The physical interpretation of the vector fields is the following. Vector field X\ indicates that l|57|l 
is invariant under time translation whereas vector fields X3 and X^ reflect the spatial invariance 
of (|57(l in x and y directions respectively. Vector field X2 demonstrates that Eq. (|57J) is invariant 
under rotation. 

The commutation relations between these vector fields become 

[X 1 ,X 2 }=0, [Vx , X 3 ] = 0, [Vx , V 4 ] = 0, [Vx , X c ] = X Ct , 

[X 2 ,X 3 ] = V 4 , [X 2 ,Xi] = -X 3 , [V 2 , X c ] = yX Cx - xX Cy , 

[V 3 ,V 4 ]=0, [X 3 ,X c ]=X Cx , [X 4 ,X c ]=X Cy , (63) 

where X Ct = c t (l - u) 2 d u , X Cm = c x (l - u) 2 d u , and X Cy = c v (l - u) 2 d u . 

Proceeding in a similar fashion for the nonintegrable case (to ^ 2) we obtain the associated 
Lie algebra as 

[X 1 ,X 2 }=0, [X 1 ,X 3 }=0, [X 1 .X 4 \=Q, 

[X 2 ,X 3 ]=X 4 , [X 2 ,X i ]=-X 3 , [X 3 ,X i ]=0, (64) 

which is obviously a subcase of ffify . 
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6 Similarity Solutions of the (2+1) Dimensional Integrable 
GF Equation 

Here we solve the characteristic equation 

dt dx _ dy du 

t £, v <t> 

and obtain the similarity reduced PDE in two independent variables. The resultant PDE is then 
analysed for its symmetry properties to reduce it to an ODE so that explicit solutions can be 
found by integrating it. In the general reduction, we obtain cylindrical function solutions. We 
derive certain physically important structures by restricting some of the symmetry parameters 
to be zero. For example, we obtain Bessel function solutions by assuming one of the symmetry 
parameters to be zero. Similarly proceeding we report certain propagating type structures and so 
on by restricting appropriately the choice of the symmetry parameters involved. In such a process, 
we report all the propagating structures discussed by Brazhnik and Tyson [1999a, b] as a subcase 
of our results. Finally, we put forth the relation between symmetries of the GF equation and the 
2-dimensional linear heat equation. 

6.1 General reductions 

Now integrating Eq. 1|65|) with the infinitesimal symmetries (|59f) along with the condition b ^ 
(as well as l,n ^ 0), we obtain the following similarity variables: 

b . 2 2\ a ■ -l ( n — bx 

Zi = - (x + y ) + ly - nx, z 2 — —t - - sin 



2 V y ; y ' b \y/P+n' 2 + 2bz 1 J ' 

u = 1 p— — . (66) 

w(z 1 ,Z2)+ J c(t,x,y)dt 

In this case also, one can obtain the linear heat equation 

Xt ~ Xxx - Xyy - X = (67) 

by assuming x — ~ (w(zi, Z2} + J c(t, x, y)dtj in i|66|) and substituing the resultant transformation 

into Eq. (|57|l . Then the similarity reduced variables l|66|l can be interpreted as giving rise to the 
linearizing transformation from the group theoretical point of view. Under the above similarity 
transformations, Eq. f^Tjl gets reduced to a PDE in two independent variables z\ and z 2 , 

2 

w Z2 + 2bw Zl + (262i + l 2 + n 2 )w ZlZl + - - Wz \ z \ +w-a = 0. (68) 

2bz\ + l z + n A 

The similarity reduced PDE (|68|) can itself be further analysed for the symmetry properties 
by once again subjecting it to the invariance analysis under the classical Lie algorithm. For this 
purpose we first redefine the variables as 

Z\ = 2bz\ + I 2 + n 2 , Z2 = Z2, w — w — a, b, n ^ (69) 

so that Eq. I|68|l becomes 

a 2 

wg 2 + Ab 2 w Zl + 4i?zxWz x z x + —Wz 2 z 2 + w = 0. (70) 

Z\ 

Eq. H70|) can be shown to possess the following symmetries 

t = c 3 , £ = 0, (fr = cxw + g(zi,z 2 ), (71) 



13 



where c\ and C3 are arbitrary constants and g{z\,z-i), which arises due to the linear nature of 

a 2 

Eq. (J7DJ, is any solution to the equation g S2 + ^b 2 g Sl + 46 2 zi<7g 1 g 1 + —gz 2 z 2 +3 = 0- The corre- 
al 

sponding Lie vector fields 



Xi = dz 2 , X 2 = wduj, X g = g{zi,z 2 )du 
lead to the infinite-dimensional Lie algebra 

[X 1 ,X 2 ]=0, [X 1 ,X g ]= X„ , [X 2 ,x g } = - 



-x 



9- 



(72) 



(73) 



where X CJz2 = g Z2 {z\ 1 z 2 )dw- By solving the characteristic equation associated with the symmetries 
(|71|l we obtain new similarity variables 



Cl£2 

£ = z\ , w — a + e V C3 



_£i _ 

/(C) + - I 9 {zuz 2 )e\ C * Z2 J d~z 2 
C3 J 



where /(£) satisfies the linear second order ODE of the form 

C 2 f" + Cf + {A + BC)f = 0, A=(a Cl /2bc 3 ) 2 , B = (1 + Cl /c 3 )/46 2 , 



(74) 



(75) 



and prime denotes differentiation with respect to C- The exact solution to Eq. (|75|l can be expressed 
in terms of cylindrical functions of the form [Murphy, 1969] 

/ = hZ&y/BC) + I 2 Z 2 (2^BC), (76) 

where Zi{C,), i — 1,2, are two linearly independent cylindrical functions and I\, I 2 are arbitrary 
constants. Then the invariant solution to the (2+l)-dimensional PDE (|57|l can be written as 



1 - a 



c 3 / (l l Z 1 {2^/BT l )+hZ 2 {2y/BT 1 )-^- J g( Zl ,z 2 )e\ c 3 Z2 J dz 2 



+ / c(t,x,y)dt 



(77) 



where z\ and z 2 are of the form (|66J) via (|69|l . Here g(z\, z 2 ) can also be taken as zero. However, 
the solution l|77|) is richer in structure when g(z~x, z 2 ) is taken as non-zero. 

6.2 Special Reductions 

Besides the above general solution, some physically interesting structures can be obtained by 
choosing some of the symmetry parameters including the arbitrary functions to be zero. In the 
following we present certain nontrivial cases. Particularly, we consider separately the two special 
cases: (a) c\ = and (b)c! = — c 3 . 



6.2(A) The choice b^O 
Case (a): c\ = 0. 

Solving the characteristic equation associated with l|71() with c\ = we obtain the similarity 
variables 



C = zi, w = a + /(C) + — / g(z 1 ,z 2 )dz 2 . 

C3 



(78) 
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Here the function /(£) satisfies the ODE 



B 



46 2 ' 



a" + f + Bf = o, 

The solution (static) to l|79|l is found to be 

/(C) = hM2y/BC) + I 2 Y (2 V / BC) 



d 



(79) 



(80) 



where Jo and Yq are zeroth order Bessel functions of first and second kinds respectively. Then the 
time-independent solution to the original PDE (|57|l becomes 



1 - a 



a + hJo h^B {b 2 {x 2 + y 2 ) + 2b(ly - nx) +P + n 2 )) 
+I 2 Y (2^/B (b 2 (x 2 + y 2 ) + 2b{ly - nx) + I 2 + n 2 )) 



(81) 



+ — / g(z 1 ,z 2 )dz2 + j c(t,x,y)dt 

The solution plot is given in Figs. (5) for g(zi,z 2 ) — and c(t,x,y) = 0. For the choice 
I\ 7^ 0, 1 2 = 0, Fig. (5a) shows a static circular symmetric pattern. If we choose I 2 ^ while 
keeping I\ — we get a static structure (Fig. (5b)) exhibiting circular symmetry with a singularity 
at the origin and it is due to the nature of Bessel function of second kind. Finally for I\,I 2 ^ 
the plot is shown in Fig. (5c). 




Figure 5: Static structures exhibiting ||%T|) circular symmetry: (a) I\ ^ 0, 1 2 — (b) Ix = 0, 1% ^ 
(c) h,I 2 ^0. 

Case (b): ci = — C3. 
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Making use of the choice c\ = 



— C3 in (|74|) we obtain the similarity variables 



C — zi, w = a + e 



/(C) + - / g(z u z 2 )e z *dz 2 
C3 



a \ 1 

2b; 



These similarity transformations reduce the PDE (|70|> to an ODE 

C 2 /" + {f + Af = 0, A 
where prime denotes differentiation with respect to £. Solving Eq. 1(83|1 we find 

f = hC mi +hC n2 , m x , 2 = ±V=A. 
Since A > 0, bounded solution can be written as 



(82) 



(83) 



(84) 



1 -a 



a + e~ Z2 h 



cos (yi log I (b 2 (x 2 + y 2 ) + 2b(ly - nx) + I 2 + n 2 ) \ ) 



sin ( VAlog I (b 2 (x 2 + y 2 ) + 2b(ly - nx) + I 2 + n 2 ) \ J 



(85) 



Here the variable z 2 is given by Eq. (|66|l . The solution plot is given in Figs. (6) for g(zi,z 2 ), 
c(t, x, y) = and it exhibits a complicated propagating pattern. 

6.2(B) The choice b = 

In the general reduction, we have assumed b 0. In this subsection, we derive certain interesting 
patterns by taking the symmetry parameter 6 = 0. For this choice the symmetries are (vide 
Eq. (29) 



r = a, £ = l, r] = n, (/) = c(t, x, y)(l - u) 2 . 
The similarity transformations 

n 

z\ = nx — ly, z 2 = ax — It, u = 1 



w(z l7 z 2 ) + J c{t,x,y)dV 



(86) 



(87) 



are found by solving the characteristic equation associated with Eq. I)86|l. Under the above set of 
transformations 187(1 , Eq. Q57J1 can be transformed to 

lw Z2 + (I 2 + n 2 )w ZlZl + a 2 w Z2Z2 + 2anw ZlZ2 + w — a = 0. (88) 

The invariance of Eq. I)88|l under the one-parameter Lie group transformations leads to the in- 
finitesimals 



n a 2 n 2 2 2 \ 

t = 2a—c\z 2 — 2— cizi + c 2 , £ = — 2a—c\Z\ + -{I +n)az 2 -c 3 , 

<j> = {c\Zi + c 4 )w + h{zi,z 2 ), 



(89) 



where w — w — a and the arbitrary function h(zi, z 2 ) satisfies the equation lh Z2 + (I 2 + n 2 )h ZlZl 



a hz 2 z 2 + 2anft, ZlZ2 +h = 0. 



Using the symmetries given in Eq. (|89|l one can construct a general similarity reduced ODE. 
However, as our aim is to show some interesting solutions we consider a lesser symmetry parameter 
group in the following. 
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(<=) 



t = " ■*''" 



f = 5 " 



t = 10 " 



* = 15 " 



t = 20 " 









Figure 6: Propagating type structures corresponding to Eq. 185JI at various instants: (a) I\,l2 7^ 
(b)7 1 ^0,7 2 = 0(c)/ 1 =Q,7 2 ^0. 



Case (a): To begin with let us consider c\ 7^ and all other parameters (02,03,04) are zero 
in 189fl . Now solving the characteristic equation associated with the symmetries we obtain the 
following similarity variables 



c 



w 



(l 2 +n 2 ) 



-z 2 — anz\Z2, 



/(C) 



Cl 



F((,z 1 )h(z 1 ,z 2 )e 



dzi 



(90) 



where F((, zi) = 
deduce an ODE 



2a 2 / 



2(1 + n )£ — a I zf . Using the above similarity transformations we 



II I 2 + n 2 
(f" + f' + Bf = 0, B=A,\1- 



2a 2 l 



Aa 2 



(91) 



where the prime stands for differentiation with respect to £• Eq. H91|) admits a time-dependent 
cylindrical wave solution of the form 



/(C) = hJ (2y/BQ + I 2 Y a (2^BQ) 



(92) 
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From l|92|) we can provide the time-dependent solution to the original PDE through the trans- 
formations (|90|l and (|87|l . The solution is plotted in Figs. (7) for various instants, exhibiting 
propagating spike like patterns. 



t = 





(4 




f = 5 " 






t = 10 




t = 15 " 



t = 20 » 







Figure 7: Bessel type propagating solutions for various instants as given by Eq. jSZ} : ( a ) 7i , 7 2 7^ 0, 
(b) 7^0,72 = 0(0) /i=0,7 2 ^0. 



Case (b): ci = 0. 

The choice ci = and c 2 , C3 , C4 ^ in Eq. (|89|l leads to the similarity transformations 



C = c 3-22 + c 2Zl, w = e^ Z2 



/(C) + - / h(z 1 ,z 2 )e kz *dz 2 

C2 



C2 



(93) 



which transforms Eq. i|88|) into an ODE of the form 

Af" + Bf + Cf = 0, A = (ac 3 ) 2 + (P + n 2 )c 2 2 + 2anc 2 c 3 , 
B = lc 3 + 2a 2 kc 3 + 2anc 4 , C = 1 + Ik + (afc) 2 , ' = 

whose general solution is given by 

/(C) = he m ^ + I 2 e m ^, m L2 - --L (b ± ^ B 2 - 4Ac) , 

and C is quadratic in time t. Substituting the similarity transformations l|93[) and l|87|) in (|95|l the 
solution to the original PDE f5?| can be obtained. 



(94) 



(95) 
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Case (c): Propagating structures. 

By assuming n — C3 = in Eq. I|94|) . the ODE becomes 



/" - hf = 0, k x = - 



(96) 



Then, the system l|57fl is found to exhibit propagating structures and the corresponding solutions 
are given by 



u = < 



1 — a < a + exp 



k (ax — It) 



h cos(v kic 2 ly) + h sin(i/ kic 2 ly) 



h(zx, z 2 ) j- Zo 



e"'~^dz 2 



+ / c(t,x,y)dt 



fci < 0; 



1 — a < a + exp 



k (ax — It) 



he VkiC2ly + he -Vkic 2 ly 



h(zi,z 2 ) ^kz? 



dz 2 



1 — a I a + exp 



k (ax — It) 



(97) 



+ / c(t,x,y)dt 



fci > 0; 



hc 2 ly + h 



h(zi,z 2 ) Js Zo 



e^dzo 



Co 



c(t, x, y)dt 



fci = 0. 



The presence of an additional function h(z±, z 2 ) in the solution (|97|l leads to a more general form. 
In particular, Eq. I|97|l exhibits the five classes of bounded travelling wave solutions reported by 
Brazhnik and Tyson [1999a] with specific choice of the parameters involved along with h(zx, z 2 ) = 
and c(t,x,y) = 0. The corresponding solutions are given below. 
The simplest travelling wave solution (Fig. 8a) 



u =1 p — -j= J: fci > 0, A = — or — 

1 + ^4 exp [fc (ax — It) ± \/kic 2 ly\ a a 



(98) 



is constructed by assuming either I\ — or I 2 = 0. When I\ = I 2 (^ 0), one obtains a V-wave 
pattern (Fig. 8b) 



1 2 
u = 1 = , fci > 0, A= -. 

1 + A exp [fc (ax — It)] cosh(\/kic 2 ly) a 



A wavefront oscillating in space (Fig. 8c) 



u = 1 



1 + A exp [fc (ax — It)] \ cos(^/kic 2 ly) \ 



A 



h 



(99) 



(100) 



exists for the choice I 2 = and fci < 0. Further for 1\ ^ and I 2 = we get a separatrix solution 
(Fig. 8d) 



u = 1 - 



1 



fci = 0, A = 



hlc 2 



1 + A\y\ exp [fc (ax — It)] ' a 

Finally, for positive fci with I\ — —I 2 the Y-wave solution (Fig. 8e) 

1 , h 

u = l — , A= — 

1 + Aexp[k (ax — lt)]\smh(^/klC 2 ly)\ a 



(101) 



(102) 



exists. In each of the above solutions A is a positive constant. As Fisher equation forms a basis for 
many nonlinear models, the above solutions are nothing but reminiscent of patterns from different 
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fields. To mention a few, V-waves are characterized in the framework of geometrical crystal growth 
related models [Schwendeman, 1996] and in excitable media [Brazhnik & Davydov, 1995], while 
space oscillating fronts are relevant to cellular flame structures and patterns in chemical reaction 
diffusion systems [Scott & Showalter, 1992; Showalter, 1995]. Besides, excitable media supporting 
space-oscillating fronts are discussed by Brazhnik et al. [1996] with a geometrical model. 




Figure 8: Five interesting classes of propagating wave patterns as obtained in ref. [Brazhnik & 
Tyson, 1999a], which follow from Eqs. (|97J) : (a) Travelling waves, (b) V-waves, (c) Oscillating 
front, (d) Separatrix solution, (e) Y- waves with £ = k(ax — It). 

6.3 Relation between symmetries of m = 2 case and of linear heat 
equation in 2- spatial dimensions 

Proceeding in a similar fashion as that for the (1+1) dimensional Fisher equation as given in 
subsection 14.31 one obtains a mapping 

zi =t, z 2 = x, Z3 = y, w = — — — (103) 

1 — u 

by solving Eqs. I|48|l for the two spatial dimensional case. Eq. (|103|) then transforms the Fisher 
equation in 2-spatial dimensions to a 2-dimensional linear heat equation 

v Xl - v Z2Z2 - v Z3Z3 - v = 0, w = l + v. (104) 

However, as in the case of the (1+1) dimensional GF equation, we have given an interpretation 
for the transformation (|103fl in terms of the Backhand transformation via Eq. Q2ip. 

7 The Nonintegrable (2+1) Dimensional GF Equation 

One can make use of the same techniques that used for the integrable case to obtain some special 
solutions for the nonintegrable case (m ^ 2). To start with, we shall consider the general reduc- 
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tions. Due to the nonintegrable nature of the equation under consideration, the general reductions 
lead to ODEs of non-Painleve type. However, we obtain certain special solutions, like plane wave 
solutions, static structures and localized structures under appropriate parametric restrictions. 



7.1 General reductions 

To begin with we derive a general similarity reduced ODE by assuming none of the parameters in 
(|60|l to be zero. The similarity variables thus obtained are 

6 2 2 o . _•■ / n — bx \ 

* = 2<* +y)+iy-™, z 2 = -*- r m- ( ^ fa+na + 2te J , 

u = w(zx, Z2). (105) 

Under the transformations H1U5|) the reduced PDE in two independent variables takes the form 

w Z2 + 2bw Zl + (2bz 1 + I 2 + n 2 )w ZlZl + 



;2 1 2\ a W Z2Z 2 1 771 



2 2 

/r,i ,72. 2\ 2 1 a W z 2 2 

(2oz\ + I +n )w Zl + 



262i + I 2 + n 2 l-w 

w + w 2 =0. (106) 



2bzi + l 2 + n 2 

Again looking the invariance properties of Eq. (|106[1 we get the following infinitesimals: 

t = ci, f = 0, = 0. (107) 
These infinitesimals reduce Eq. (|106J) to an ODE 

46 2 ( U" + f + ^-J' 2 ] +f-f = (108) 



i-r 

under a set of similarity transformations 

C = zi, w = f((), zi = 2bzi + I 2 + n 2 . (109) 

As Eq. I|1U8|I is not integrable, we look for special reductions with lesser number of parameters 
involved. 

7.2 Special reductions 

For the special case 6 = the infinitesimals obtained are r = a,£ = I, r\ = n, <fi = 0. The 
associated similarity transformations Z\ — nx — ly, z-i — ax — It, u — w(zi, Z2) reduces the Fisher 
type equation (|57|l to a PDE 

lw Z2 + (I 2 + n 2 )w ZlZl + a 2 w Z2Z2 + 2anw ZlZ2 + 



1 — w 

x[(l 2 +n 2 )w 2 Zl +a 2 w 2 Z2 +2anw Zl w Z2 ]+w-w 2 =0 (110) 

in two independent variables z\ and Zi- Applying classical Lie algorithm on Eq. (|110|l one obtains 
the following infinitesimals 

njfii +c 2 ) 

T = C 2 , , = 0. Ill 

a 

The corresponding similarity variables, £ = — ac2Z\ + n(c\ + C2)z2 and w — /(C), reduce Eq. H110JI 
to 

Af" + Bf ^f 2 - f + f 2 = 0, (112) 

with / = 1 — /, A = a 2 (c 2 n 2 + c 2 Z 2 ), B = nl{c\ + c 2 ) and prime stands for differentiation with 
respect to C- Because of the nonintegrable nature of Eq. (|112f) . we look for subcases by assuming- 
one or more of the vector fields to be zero. In the following we report some of the nontrivial cases. 
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7.2(A) Plane wave solutions 

For m = 0, making use of a similar procedure as in the previous case we can obtain the plane 
wave solution. That is, m — leads to the ODE for £ = C/VA 



VA 



(113) 



It is then straightforward to check that for Bj\J~A — 5/V6, Eq. 1(113(1 satisfies the Painleve property 
so that the solution to the original PDE is found to be 



u(x,y,t) = 1 - 



k I n(ci + C2)(ax — It) — acMy — nx) 

1 + vl cxp 7u 



(114) 



Here we obtain a propagating plane wave and the pattern is plotted in (Figs. 9). 

m m ® m 




Figure 9: An exact propagating plane wave solution Eq. (|114J) at different instants of time: (a) 
t = 0, (b) t = 5, (c) t = 10, (d) t = 15. 



7.2(B) Static and localized structures 

Substituting c\ = — c 2 in (|lll|l . we obtain the similarity transformations £ = Z\, w — /(C)- This 
transforms Eq. I|11U|) to an ODE 



/" - jf' 2 ~ Af(l - /) = 0, A 



1 



I 2 + n 2 ' 



(115) 



For Eq. (|115fl . in addition to the elliptic function solutions of the form tabulated in Table 1, we 
obtain a particular planar solitary wave solution (Fig. 10) 



u = 1 



3 - 2m 
2 -2m 



sec h 2 I 2 - 



.4 



1 — TO 



, TO < 1. 



(116) 



We note the solutions given in (|114(l and (|116J) are 2-dimensional generalizations of 1)36(1 and 
(l4*21l respectively. 

Besides the above, the choice 6 = 0, a — 0, I — reduces the (2+1) dimensional Fisher 
equation 157(1 to that in (1+1) dimensions for to ^ 2. 



8 Nonclassical Reductions 

Generally, for many PDEs there exist symmetry reductions that are not obtained by using the 
classical Lie group method. As a consequence, there have been several generalizations of the 
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Figure 10: Fig. 10. A planar solitary wave solution for Eq. (|115fl in the particular limit. 



Table 2: Elliptic Function Solutions of Eq. QgJ for 



1/2. 



S.No. 



Order of Roots 



Function v(z) 



modulus k 



ci > c 2 > > v 



Cl > c 2 > > v 



ci > v > c 2 > 



ci > v > c 2 > 



V = Cl 



Cl 



sn2 -^-( z -^)> k 



< ->.sn 2 [ ^~{z - h),k 



ClC 2 



ci - (ci - c 2 )sn 2 
v = ci - (ci - c 2 )sn 2 




, 2 _ Cl - C 2 
Cl 



, 2 = Cl - C 2 
Cl 

, 9 Cl - C 2 



fc 2 = 



Cl 

ci - c 2 



Cl 



classical Lie group method for symmetry reductions. Bluman and Cole [1969], then proposed the 
so-called non-classical method of group-invariant solutions in the study of symmetry reductions 
of a linear heat equation. An algorithm for calculating the determining equations associated with 
the nonclassical method was presented by Clarkson and Mansfield [1993]. This procedure has 
been applied to several nonlinear systems and in some cases, such as the Boussinesq equation, the 
Burgers' equation and the FitzHugh-Nagumo equation, new similarity reductions not obtainable 
through classical symmetries have been found [Levi & Winternitz, 1989; Arrigo et al. 1993; 
Nucci & Clarkson, 1992] (see also Mansfield et al. [1998] and references therein). In the classical 
method we consider the infinitesimal generator, associated with the one-parameter Lie group of 
transformations, 



V = £(t, x, u)d x + r(t, x, u)dt + 4>(t, x, u)d u 



(117) 



which leaves the system <|13t invariant. But in the nonclassical method, one requires the given 
equation l|13|l and the surface condition 



£(t, x, u)u x + rit, x, u)u t — 4>(t, x, u) = 



(118) 



together to be invariant under the transformation with the infinitesimal generator (|117|l . In this 
case one may obtain a larger set of solutions than that for the classical method. Moreover, 
significant progress has been made in the study of nonclassical symmetries for nonlinear PDEs of 
diffusive type. To mention a few, Gandarias and Bruzon, [1999] have obtained several solutions 
for a family of Cahn-Hilliard equations that are not invariant under any Lie group admitted by 
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the equation. In a similar fashion, separable new solutions that are not obtainable via classical 
method have been reported for a mathematical model [Gandarias, 2001] of fast diffusion through 
nonclassical method. This prompts us to look for nonclassical symmetries associated with the GF 
equation ijjjj. 



8.1 The (1+1) dimensional Fisher equation 

There are two types of nonclassical symmetries: those where the infinitesimal r is non-zero and 
those where it is zero. In the first case, we can assume without loss of generality that r = 1, 
while in the second case we can assume £ ^ and r = 0. In the following we investigate them 
separately. 



8.1(A) The case t^O 

We set t — 1 in the invariant surface condition without loss of generality and use it together with 
its differential consequences to eliminate u t and so Eq. I|13[) takes the form 

u xx = 4> - £u x - — u 2 x - u + u 2 . (119) 
1 — u 

Applying the classical Lie algorithm to l|119|) and eliminating the highest derivative involving the 
variable x, the coefficients of the linearly independent expressions in the remaining derivatives are 
set equal to zero. From the resultant set of determining equations for £, t, and 0, which are in 
general nonlinear, we get 

r=l, £ = cj) = d{t,x)(l-u) 2 , m = 2, (120a) 
r=l, £ = b, = 0, m^2, (120b) 

where b is an arbitrary constant. In Eq. l|12Uall . d(t, x) satisfies d\ ~ d xx — d = 0. 
The corresponding similarity variables for the m = 2 case are 

z = x-bt, u=l — ^— — , (121) 

w(z) + J d(t,x)dt 

which are similar to the ones obtained by the classical method earlier in Sec0] Thus we are lead 
to the similarity reduced differential equation 

w" + bw' + w = 0, w = W-l, ' = d/dz, (122) 

which is similar to (|31|l and the results then follows as before. Similar is the case with m ^ 2. 
Thus both classical and nonclassical reductions lead to the same solutions for the present system. 



8.1(B) The case r = 0, £ ^ 

Let t = and £ ^ 0. We then set £ = 1, without loss of generality. The invariance surface 
condition simplifies the Fisher equation to 

m 

u t - <t> x - 4>4>u - z (j) 2 ~u + u 2 = 0. (123) 

1 — u 

Applying the classical Lie algorithm to the above equation we end up with a more complicated 
nonlinear PDE 

4>t - 4>xx - ^4>4>xu - 4> 2 4>au + u(i - u)4> u + (2u - \)4> 

r0 3 =0. (124) 



1 — u 1 — u (1 — u) 2 



24 



This system is considerably more complex than the system (lKi[l and hence cannot be solved in 
general. We have made different ansatze for 4>(t, x, u) and found that all of them lead to only the 
group invariant solutions found by the classical algorithm. 

However, for the m = case, by assuming <fi to be independent of x and t we get 4> — 

mW- 1, t = 0, £ = 1 which is different from that obtained through the classical method. 

However this choice leads to a singular solution of the form 




(125) 



8.2 The (2+1) dimensional GF equation 

We extend the nonclassical symmetry analysis to the (2+1) dimensional Fisher equation (|57|l and 
in this analysis there are three cases to consider, namely, (i) t / 0, (ii) r = 0, rj ^ 0, and (hi) 

T = o, t) = o, a- 0. 



8.2(A) The case r ^ 

We use the invariance surface condition 



u t = (j>- £u x - i]u y (126) 

with its differential consequences to eliminate u±. Applying the classical Lie algorithm to the 
Fisher equation 

4> - (,U X - TjUy - U XX ~ Uyy ~ —— (v^. + Uy ) - u + = (127) 

and then equating coefficients of various powers of derivatives of the dependent variable u we have 
a set of determining equations. Solving them we get finally 

t = 1, £ = a^y + <24, -q — — a^x + C4, (f) — d(t, x, y)(l — it) 2 , m = 2, (128a) 
t = 1, £ = a^y + 04, 7/ = — 032; + C4, = 0, m 7^ 2. (128b) 

Here d(t, sc, y) satisfies the equation <f t — — d yy — d = and 03, 04, C4 are arbitrary constants. 
Comparing this set of symmetries with that obtained by the classical method we find that they 
are similar to the latter case. So the similarity reductions will reduce the PDE to the same ODE 
as that obtained by the classical method. 



8.2(B) The case t = 0, rj ^ 

Here we set rj — 1, without loss of generality, and analogous to the above procedure we use 
the above choice to eliminate all y-derivatives. Then applying the classical Lie algorithm to the 
resulting equation we obtain a more complicated nonlinear PDE. Solving the resultant system of 
determining equations we obtain the following symmetries 

(*) r = 0, C = ci, 17 = 1, = 0, (129a) 

(ii) r = 0, £ = ±4|, = 0, (129b) 

for all values of m. The similarity reduction variables z\ — c\y — x, Z2 = t, u = w{z\, Z2) associ- 
ated with (|129af> reduces the GF equation in 2-spatial dimensions to that in 1-spatial dimension. 
Similarly, the similarity transformations Z\ — ±V2y — x, z% = t, u = w (zi,z 2 ) for IT29bT) again 
reduces the original PDE to that in 1-dimension and hence the results follow as in the case of the 
classical method. 
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8.2(C) The case t = 0, r) = 0, £^0 



In this case we set £ — 1, without loss of generality, and as above we make use of the invariance 
condition u x = <f> to eliminate all ^-derivatives. Then we apply the classical Lie algorithm to the 
resulting equation and we get a nonlinear PDE of the form 



(l-u)L(l-u) 



+ Uy) + (1 - 2u + 24> ux )4> + 2(j) uy u y - u(l - u)(j> u 

{ul + (j) 2 ) + <j> u {ul + 4> 2 ) + 2(f>yUy + 2# x 



0. (130) 



Solving the system of determining equations we obtain <p — V m. Thus the (2+l)-dimensional 
GF type equation, under the similarity transformation z% = y, Z2 = t, u = w(z\,Z2) reduces to 
that in (l+l)-dimensions and hence the results as in the case of the classical Lie algorithm case 
follow. 



9 Generalized Conditional Symmetry Reductions 

Recently Fokas and Liu [1994a, b] proposed the notion of Generalized Conditional Symmetry (GCS) 
and applied it to construct some physically interesting exact solutions of certain nonlinear non- 
integrable PDEs. Such exact solutions are of primary importance because they identify certain 
interesting and novel physical phenomena and moreover such solutions may be hard to identify 
or may not be quite transparent from the numerical solution of a nonlinear PDE. In this method 
the PDE is reduced to an ODE in terms of certain types of generalized conditional symmetries. 
Indeed the GCS can be considered as a natural generalization of the nonclassical method just as 
the generalized symmetry method [Olver, 1986] is considered as a generalization of the Lie-point 
symmetry method. Particularly, for the reaction diffusion type equations, several new separable 
solutions that cannot be obtained via the non-classical and non-local symmetry methods have been 
derived through this method [Qu, 1997; Qu, 1999a,b; Chou & Ku, 1999; Qu et a/., 2000]. Further, 
it is recently shown that the GCS method is the most efficient tool for solving the problem of 
dimensional reduction of initial value problems for evolutionary type PDEs in a purely algebraic 
way [Zhadanov & Andreitsev, 2000]. Motivated by these facts we investigate the existence of 
generalized conditional symmetries, if any, for the GF equation. 

9.1 The method 

Let K(t,u) denotes a function that depends on a differentiable manner of t,u,u x ,u xx , . . .. The 
function a(t, x, u) is a generalized symmetry of the evolution equation u t = K(t, u) iff 

^ + [K,a]=0, (131) 

when u t — K(t,u), where [K, a] = a'K — K'a, and the prime denotes the Gateaux derivative. 
The function a(t, x, u) is a GCS of equation (|13|) if there exists a function A such that 

-£- + [K,a]=A(t,x,u,o), A(t,x,u,0) = 0, (132) 
ot 

where K(t, u), and a(t, x, u) are differentiable functions of t, x, u, u x , u xx , . . . and A(t, x, u, a) is a 
differentiable function of t, x, u, u x , u xx ,..., and a, a x , a xx , ... It is obvious from l|132|) that Eq. (|131f) 
admits a GCS a iff 

a'K | CT=0 = 0, (133) 
provided a is explicitly independent of time t. 
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9.2 The (1+1)- dimensional GF equation 

Let the Fisher equation admits a GCS of the form 

a = u xx + H(u)u x + F(u)u x + G(u), (134) 
where H(u), F(u) and G(u) are arbitrary functions u. Substituting l|134|) into (|133[) one gets 



Huu — 2H(2H U — H 2 ) 



{l-uf 



1 -u 



- m 



l-u 



(H u - 2H 2 ) 



F U u — 2HF U — AF(H U — H 2 



-(F U +4FH) 



AmF 



l-u"' - ' ' (l-u) 2 

2F{F U - HF) - AG(H U - H 2 ) + (2u - l)H - u(l - u)H u 
5mG 



1 - u 



(AHG - 2F 2 + G u ) 



{l-uf 



(u(l - u) + 2G) F u 



AFG 



l-u 



-H 



u x - u(l - u)G u + (I - 2u)G 



-2G 2 (H 



-) = 0. (135) 



l-u' 

Equating coefficients of various powers of u x to zero and then solving the resulting overdetermined 
equations we get 

m 



H(u) 



G(u)=u(l — u), F(u) = ci, V m, 



l-u' 

where ci and C2 are arbitrary constants. The corresponding GCS are 



a = u a 



-U x + C\U X + u — u . 



(136) 



(137) 



l-u 

With the above some exact solutions can be constructed. The ODE a — is solved to obtain u 
as a function of x with integration constants dependent on time t alone. Then substituting this 
solution into the governing equation, the time evolution of these time-independent constants are 
determined. The exact solutions are then obtained by solving these systems. 
The GCS (H32J leads to 

to 



1 - u 



u' 2 + ciu' 



u-u A = 



which cannot be integrated as such except for to = 2. Thus for m = 2, we obtain 

-ci ±y/ct-4 



1 - [1 + A(t)e miX + B(t)e r 



mi, 2 



Substituting (|139|) into the Fisher equation (|13fl we get 



mi o 



Ci ± yjc\ - 4c 2 



(138) 



(139) 



(140a) 



(140b) 



for to = 2. We wish to mention that the solution l|140a(l can be recovered from the linear heat 
equation Wang et al. [1996] and the results are similar to that obtained through the classical 
algorithm. Further Eq. (1 140 b() is similar to Eq. I|33|l and hence we obtain the various propagating 
structures already discussed through the classical method. 

Eq. I|138|l is of non-Painleve type for all other values of m 2) and is similar to (|35|l and 
hence the previous results follow. Thus GCS reductions lead to the same results as that obtained 
through the classical method for the modified perturbed GF equation. 
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10 Conclusions 



In this paper, we have investigated the intcgrability/symmetry properties of the generalized Fisher 
type nonlinear reaction-diffusion equation in both (1+1) and (2+1) dimensions. The singularity 
structure analysis singles out the m = 2 case as the only system parameter for which the GF 
equation is free from movable critical singular manifolds. Further a Backlund transformation for 
the intcgrable case is shown to linearize the GF equation to the linear heat equation. The symmetry 
analysis shows that the system under consideration possesses interesting Lie point symmetries 
that lead to the infinite-dimensional Lie algebra for the integrable case, thereby exhibiting various 
interesting patterns and dynamics, and giving rise to a group theoretical interpretation for the 
system. In addition to the above, we have given exact and numerical solutions for specific choices of 
the nonintegrable case. Our studies on the GF equation reveals the existence of a large number of 
interesting wave patterns, static and localized structures. We have also carried out the nonclassical 
and generalized conditional symmetry reductions. 
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